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Summary 


Recent  work  on  the  short-term  statistics  of  scintillation  is 
extended  to  a more  detailed  investigation  of  the  effective  number  of 
independent  turbulence  regions  on  a propagation  path  as  a function  of 
wind  direction  and  other  parameters  The  resultant  effect  on  the  spread 
in  strength-of-f luctuat ion  measurements  vs.  averaging  time  is  discussed, 
including  both  microtliermal  and  scintillation  data  and  their  inter- 
relationship. These  considerations  are  related  to  spatial  correlation 
measurements  of  the  turbulence  strength  (mlcrothermal  envelope)  as  a 
function  of  averaging  time,  including  the  Influence  of  pronounced  inter- 
mittency  of  turbulence. 

In  a related  discussion,  the  conditional  statistics  of  scintillation, 
based  on  an  initial  comparison  of  the  instantaneous  target  irradiance 
with  a predetermined  threshold  level,  are  investigated  as  a function  of 
elapsed  time.  With  tlie  assumption  that  the  two-point  log  scintillation 
process  is  a bivariate  normal,  the  problem  is  readily  solvable  using 
empirically  derived  correlation  functions.  The  validity  and  limitations 
of  .his  assumption  are  briefly  considered,  and  representative  results 
are  given. 

New  multiwavelength  facilities  are  described  which  will  be  used  for 
the  completion  of  the  experiments  on  mean  target  irradiance  and  its 
fluctuations,  with  and  witliout  cancellation  of  turbulence- induced  beam 
wander.  Other  near-term  plans  are  also  briefly  reviewed. 
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I.  INTRODUCTION 

This  report  will  briefly  summarize  recent  progress  and  ongoing 
efforts  on  several  topics  in  optical  propagation  through  atmospheric 
turbulence.  In  Section  II,  ue  describe  some  further  results  on  the 
slio»"t-term  statistics  of  turbulence  and  scintillation,  which  supple- 
ments tlie  detailed  treatment  given  in  the  preceding  report  on  this 
program.  In  Section  III,  we  discuss  a new  topic,  the  conditional  fading 
statistics  of  a scintillating  signal  or  target  illumination,  in  which 
the  state  of  the  "channel"  is  investigated  immediately  following  the 
onset  of  a "good"  interval.  Finally,  in  Section  IV  we  describe  new 
1 experimental  work  and  facilities,  including  mult iwave] ength  target- 

irradlance  investigations  with  beam-wander  cancellation,  very-long- 
path  or  strong-scattering  scintillation  phenomena,  and  other  areas  of 
Interest . 

I 

II.  FURTHER  RESULTS  ON  SHORT-TERM  STATISTICS 

In  the  preceding  technical  report  on  this  program,^  we  presented 
a complete  theoretical  treatment  of  the  short-term  statistics  of  tur- 
bulence and  scintillation,  including  data-spread  and  confidence- 
interval  considercitions  as  a function  of  averaging  time.  In  this  sec- 
tion, we  supplement  the  previous  discussion  with  representative  results 
of  recent  experiments  on  this  topic.  Our  related  work  on  computer 
simulation  of  the  short-term  propagation  problem  will  be  described  in  a 
subsequent  report. 

A.  Relationship  Between  Microthermal  and  Scintillation  Data  Spread 

In  a discretized  model  of  the  propagation  path,^  we  idealize 
the  turbulence  into  N elements,  the  strength  of  turbulence  being 
totally  correlated  throughout  one  element  and  totally  uncorrelated  between 
any  two  elements.  The  length  of  an  element  depends  on  the  averaging  time 


1.  J.  R.  Kerr,  et  al,  "Propagation  of  Multiwavelength  Laser  Radiation 
Through  Atmospheric  Turbulence,"  RADC-TR-74-320,  November,  1974. 
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n 


variance  of  n 


(t)  used  in  detcrmininn  "stren,;/th  of  turbulence"  / C • and  on  the 

V "J 

direction  of  the  wind  relative  to  tlie  path.  The  variance  of  a 

sot  of  measurements  of  the  scintillation  lo^  amplitude  variance  ^Var 

is  then  related  to  the  variance  of  mi c rot herma 1 measurements  ^Var  \ 

\ ) 


Va  r a 2 

T 


Var  C‘ 


The  quantity  f i*’  derived  from  N,  the  number  of  correlation 

lengths  of  along  the  path,  and  the  spatial  weighting  function  (f) 

determined  ^ by  beam  geometry.  If  the  log  amplitude  variance  is 
written  ns 


f . C ^ 
1 n , 


then  we  have 


(s  q) 

N ,,  \i=l  / 

eff  = — 

E 


< N 


For  a spherical  wave,  f is  a discrete  approxim.at ion  to 


5/6  , 

f(z)  a.  (f-) 
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where  z is  the  path  variable  and  L is  the  total  pathlength.  In  the 

case  of  e.g.  focused  beams  or  other  systems  yielding  more  localized 

weighting  functions,  N becomes  much  less  than  N and  Var 

eft  X 

increases  accordingly.  In  the  limit  of  f as  a delta  function,  the 

variance  in  o ^ for  any  given  t will  be  the  same  as  that  for  . 

f n 

For  averaging  times  greater  than  the  integral  scale  of  die  ^ 

squared  microtherm;,il  fluctuations  (l._2),  the  normalized  variance  of 

,2  ' * 

measurements  will  be  inversely  proportional  to  t (Ref.  1).  The 

corresponding  dependence  for  the  normalized  variance  of  o ^ depends  on 

T 

the  wind  field.  If  we  extend' our  idealization  to  assume  Jarge-scale 
regular  flow  with  "frozen-in"  strength  of  turbulence  over  scales 
comparable  with  the  pathlength,  we  may  distingiKsh  the  following  two 
cases : 

1.  When  the  wind  is  parallel  to  the  path,  the  z-correlat ion 

length  of  is  proportional  to  t,  which  cancels  the 

usual  effect^of  increased  averaging,  and  thus  Var  o ^ is 

_ I T 

indeoendent  of  i.  That  is,  N t , to  be  cons’  .ent  with 
F.q.  1. 

2.  When  the  wind  is  perpendicular  to  the  path,  the  z-cor- 

, _ 1 

relation  length  is  Independent  of  t,  Var  t , and 

N is  independent  of  t. 

When  we  compare  this  useful  but  idealized  conceptual  model  to 

experimental  data,  subtleties  and  complications  immediately  arise,  as 

shown  by  the  following  examples.  In  Table  I,  we  show  the  basic 

parameters  of  nine  .experimental  runs  comparing  normalized  and 

n 

variances  vs.  averagi\g  time.  Runs  1-6  involve  winds  basically  parallel 
to  the  path,  and  the  results  for  (Eq.  1)  are  shown  in  Fig.  la. 

According  to  the  foregoing  reasoning,  we  would  ideally  expect  a slope 
of  +1  for  these  curves,  and  a tendency  to  break  in  this  direction  may 
be  noted.  In  Fig.  lb,  curves  are  given  for  wind  directions  partially 
or  entirely  perpendicular  to  the  path;  two  of  these  tend  to  be  flat, 
as  expected.  Taking  into  account  all  conditions  and  averaging  times, 
values  of  were  seen  to  range  from  approximately  8 to  200. 
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TABLE  T 


Wavelength:  A880X 

Pachlength:  1.6  km 

Beam  Geometry:  Approximate  Spherical  Wave 


Parameters 

for  Data  on 

Var  C 2 

and  Var  0“^  Shown 

T 

in  Figs.  1,  2, 

Wind  Direction 
Relative  to 

1 

Run  No . 

Symbol  on 

Figs. 

V (m/sec) 

Optical  Path 

0^ 

1 

0 

3.5 

0° 

0.31 

2 

□ 

<1 

0° 

0.24 

3 

0 

3 ; 

0° 

0.28 

4 

A 

3.5 

180° 

0.30 

5 

+ 

9 

-150 

• 0.076 

6 

3.5 

195" 

i 

-45° 

0.31 

7 

• 

4 

0.25 

8 

A 

0.5-2 

90° 

0.33 

9 

■ 

1.5 

90° 

« 

0.32 

-6- 


T (Si'c) 


■fne  reciproial 
ri'f;ions  / J 

V 


niinibi  r 


) 


aloiiK 


of  ef  f ect  i ve ly-indepc'ndont  tiirbulonco 
a propagation  path,  vs.  aviTasinp  tinu 


(').  Ihf  runs  aro  identified  by  tlie  symbols  in  Table  I. 
Wind  parallel  to  path. 


1 


5 


10 


20 


60 


Figure'  lb 


i 


1 (see) 


The  reciprocal  number  of  effect ively-independent  turbulence 
regions  / J_\  along  a propagation  path,  vs.  averaging 


time  (i).  The  runs  are  identified  by  the  sumbols  in  Table  1. 
Wind  components  normal  to  the  path. 
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We  now  examine  the  results  for  Var  alone  (Fig.  2a).  The  slopes 
are  in  general  less  than  the  expected  ^ value  of  (minus)  one,  suggest- 
ing that  T is  not  in  general  >>  If  the  wind  speed  is  v,  the 

value  of  vl^^2  is  a function  of  the  degree  of  spatial  intermittency  of 
the  turbulence  (Ref.  1),  and  in  some  of  the  curves  may  be  indicated  by 
a breakpoint  in  Fig.  2a. 

The  results  for  Var  likewise  do  not  show  the  simple  behavior  of 
the  idealized  model  (Fig.  2b).  We  point  out  three  complicating  factors: 
1.  A real  wind  field  may  be  expected  to  be  irregular  over  the 
scales  of  Interest  here,  so  that  regions  of  given  turbulence 
strength  are  not  "frozen-in".  This  will  cause  a tendency 
for  a (-1)  rather  than  (0)  slope  for  parallel  wind  in  Fig.  2b. 

2-  If  1 is  not  >>  I^j2>  if  fhere  are  large-siale  correlations 

!for  low  spatial  frequencies  in  AT“^ , the  slopes  in  Fig.  2b  will 
tend  towards  zero.  This  will  be  especially  noticeable  for 
transverse  wind,  where  they  would  otherwise  be  (-1),  and  this 
is  evident  for  runs  7 and  9. 

3.  For  small  t the  scale  of  the  actual  scintillations 

themselves,  Eq.  14a  of  Ref.  1 will  apply,  and  slopes  will  tend 

to  (-1)  as  in  Fig.  2b.  However,  these  shorter  averaging  times 

are  not  generally  useful  for  indications  of  C ^ or  o*^. 

n 

B . Microthermal  Envelope  Correlations 

These  subtleties  can  be  further  explored  using  spaced  micro- 

! thermal-probe-pair  measurements  of  correlations  vs.  averaging  time. 

The  measurement  scheme  is  described  in  Fig.  13  of  Ref.  1.  In  Figs.  3a, 

I b correlations  of  the  mean-square  fluctuations  AT^  | (-v  C are  given 

^ for  three  values  of  t,  for  separations  nearly  parallel  to  the  wind 

direction  (the  line  along  which  the  respective  probe  pairs  are  situated 
is  slightly  offset  from  the  wind  vector  to  preclude  wake  effects  from 
the  instruments).  The  basic  spatial  integral  scale  of  AT^  can  be  deduced 
from  the  lowest  or  t = 0.4  sec  curve,  and  the  tendency  of  the  integral 
scale  to  become  equal  to  v t/2  for  larger  t is  evident  from  the  other 
two  curves. 


O 


I 

< 


i 
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Similar  curves  for  a transverse  wind  are  given  in  Figs.  4a,  b.  The 
tendency  for  the  large-T  curve  to  decorrelate  at  larger  separations  is 
evident  in  Fig.  4b.  In  Fig.  4a,  an  anomaly  is  observed  whereby  the  basic 
squared-mlcrothermal-f luctuations  exhibit  a correlation  tail  at  large 
separations,  which  is  also  manifested  in  the  large-i  case. 

C.  A Quantitative  View  of  Turbulence  Intermittency 

A simple  operational  measure  of  the  "intermittency"  of  turbulence 
is 

Standard  Deviation  for  C ^ 

. = ^ 

C 2 
n 

where  t may  be  taken  as  e.g.  10  seconds.  In  Figure  5 we  show  differ- 
ential probe  measurements  of  AT  for  values  of  R of  0.46  and  0.68,  which 
might  be  described  as  intermediate  and  pronounced  intermittency,  respec- 
tively. 

In  Table  II,  data  ranging  from  cases  of  light  intermittency  (R  = 

0.26)  to  heavy  intermittency  (R  = 0.77)  are  shown.  The  integral  scale 
of  AT2  after  filtering  or  smoothing  with  t = 1 sec,  is  the  integral 
of  the  autocorrelation  function  of  this  smoothed  process;  it  is  seen  to 
increase  with  increased  intermittency,  as  expected,  and  decreases  toward 
t/2  only  for  light  intermittency.  This  indicates  that  intermittency 
does  indeed  introduce  a long  tail  in  the  correlation  function  of  tlie 
basic,  unfiltered  (AT^)  process.  Similarly,  the  flatness  factor  62 

A T 

(Ref.  1)  associated  with  the  probability  distribution  of  the  un- 
filtered fluctuations  is  seen  to  be  much  higher  in  the  intermittent 
case,  directly  indicating  the  predominance  of  occasional  excursions  of 
extreme  magnitude  in  that  case.  It  may  be  also  noted  that  intermittency 
is  generally  associated  with  low-wind,  poorly-mixed  conditions  charac- 
terized by  occasional  "bursts"  of  turbulence. 


UiJ. 


III.  CONDITIONAL  FADING  STATISTICS 


A topic  related  to  the  short-term  considerations  discussed  pre- 
viously is  that  of  the  conditional  statistics  of  scintillation  or  irra- 
diance  fluctuations.  The  applications  of  this  concept  include  "pilot- 
tone"  or  burst  communications  systems  in  which  the  channel  is  monitored 
and  used  only  when  it  is  "good",  and  target  illumination  systems  in 
which  e.g.  a high-power  laser  is  utilized  immediately  after  a low-power 
laser  (or  some  other  channel  monitoring  scheme)  indicates  a favorable 
short-term  path.  The  applicable  question  may  be  posed  as  follows: 
given  that  the  instantaneous  target  irradlance  is  above  some  predeter- 
mined threshold,  what  is  the  probability  distribution  of  the  irradlance 
at  an  arbitrary  time  t later?  The  present  discussion  is  limited  to 
the  linear  propagation  case. 

A.  Analysis  of  the  Conditional  Problem 

In  previous  work  on  this  program,  it  has  been  determined  that 
scintillation  does  not  constitute  a Markov  process.  It  has  been  deter- 
mined by  ourselves  and  a number  of  other  investigators,  however,  that 
this  phenomenon  is  for  most  practical  purposes  log-normally-distributed. 

To  the  extent  that  we  can  assume  that  the  two-point  log  distribution  is 

2 

bivariate  normal,  the  problem  is  readily  tractable.  However,  we  point 

out  that  the  conditional  statistics  can  be  very  sensitive  to  departures 

from  this  distribution  and  theoretical  arguments  suggest  that  exact  log 

3 

normality  is  not  even  possible  for  the  one-point  distribution.  Fur- 
thermore, as  described  below,  our  preliminary  analysis  of  experimental 
data  shows  significanc  departure  of  the  two-point  distribution  from  the 
normal  under  conditions  of  strong  scintillation. 

Nevertheless,  the  basic  approach  given  here  to  the  conditional 
problem  can  be  numerically  extended  to  utilize  any  (non-gaussian)  joint 


2.  S.  0.  Rice,  "Mathematical  Analysis  of  Random  Noise",  in  Selected 
Papers  on  Noise  and  Stochastic  Processes,  N.  Wax.,  ed.,  Dover  Pub- 
lications, Inc.,  N.Y. , 1954. 

3.  J.  W.  Strohbehn,  T.  Wang,  J.  P.  Speck,  Radio  Science  j_0,  59-70, 
January  1975. 
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probability  distribution  which  may  be  empirically  or  otherwise  deter- 
mined. We  will  confine  this  discussion  to  the  log-normal  case  for 
simplicity,  and  empirical  verification  under  certain  conditions  will 
be  given  below. 

Let  us  therefore  assume  that  the  two-point  or  joint  probability 
distribution  of  the  log  irradiance  is  bivariate  normal.  The  important 
parameter  is  the  correlation  coefficient  p(t),  which  is  a function  of 
the  time  delay  t separating  the  initial  and  final  observations  of  the 
irradiance.  The  joint  distribution  of  the  log  irradiance  x is 


Pr(x(t)  ,x(t  + t))  = 


exp  < - 


1. 

,2  / 1 _.2 I 


2o‘^(1-P‘(t)) 

o 


.^X^(t)  - 2p(T)x(t)x(t  + t)  + x^(t  + t))  I 


(6) 


where  x i®  normally  distributed  with  variance  and  zero  mean, 

and 

p(,)  , 

0^ 

o 


The  conditional  distribution  of  x(t  + ^),  given  that  x(^i  = L,  is 


Pr(x(t  + t) I x(t)  = L)  = 


(l-p'^(t)) 
o 


1/2 


exp  / 


-1 


202(1-P^(t)) 

o 


• (x(t  + t)  - p(i)L) 


2 '> 


(7) 
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which  is  normal  with  the  following  mean  and  variance; 


Mean  = Lp(t) 

Var  = 0^  (l-p2(T))  . (g) 

We  are  now  interested  in  the  probability  that  x(t  + t)  is  less 

than  a fixed  "fade"  level  (^i),  given  that  x(t)  exceeds  a threshold 

"pilot  level"  (5,  ) : 
o 


Pr(x(t  + t)  < I x(t)  > l^)  =J  dx(t)j*  d(x(t  + t))  P(x(t),x(t  + T)) 


X 

I 


dX 


J 2o2 

o/ia 


(9) 


We  define  the  following  quantities:^* 


1 -x-^ 

Z(x)  = e 

/2^ 


(9a) 


/' 


P(x)  = I Z(t)dt 


(9b) 


00 

/ 


Q(x)  = I Z(t)dt 
'x 


(9c) 


4.  M.  Abramowitz  and  I.  SLegun,  ed..  Handbook  of  Mathematical  Functions. 
Dover  Publications,  Inc.,  N.Y.  931-940,  1965. 


g(x,y,p) 


(2tt/w2)-1  j"  _ (x2-2pxy  + y2) 


(9d) 


L(h,k,p) 


00  00 


■J  “■‘i 


dy  g(x,y,p) 


(9e) 


We  note  that 


P(x)  + Q(x)  = 1 
P(-x)  = Q(x) 


00  k 

L(h,-k,-p)  = J dxjdy  g(x,y,p) 
h 


With  these  definitions,  we  have 


lO 


(10a) 

(10b) 

(10c) 


(11) 


This  expression  must  be  evaluated  numerically  for  various  values  of  £ 

0 

£l,  a , and  p.  The  calculation  is  expedited  by  use  of  the  formula 
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L(h,k,p)  = Q(h)Q(k)  + 


Z^^^(h)Z^^^(k)  n 

(tl  + 1)  ! ^ 


+ 1 


(12) 


n=o 


where 


^ Z(x)  = -xZ^'^  “ ^^(x)  - (n  - 1)  Z^"  “ ^^x)  (13) 

dx 


In  practice,  in  the  worst  cases  about  50  terms  must  be  retained  in  the 
sum  for  the  accuracy  to  be  better  than  1 x 10 

We  have  performed  computer  runs  for  several  different  pilot  and 
fade  levels,  using  correlation  functions  (Fig.  6,6a)  obtained  empirically 
under  three  different  scintillation  conditions  at  10. 6p: 

Strength  of  Scintillation  Wind  Speed  Symbol 


Run 

A 

a Va  = 0-307 

1 m/s  variable 

A 

Run 

B 

0 Va  = 0.010 

3.5  m/s  at  A5°  to  path 

0 

Run 

C 

a 7a  = 0.0015 
0 ' ^ 

9 m/s  at  A5°  to  path 

X 

The  pilot  and  fade  levels  are  conveniently  normalized  by  the  mean  irra- 
diance  I.  We  thus  calculate  the  probability  that  In  ^I(t)/iJ  < ln(I^/l) 
given  that  In  [l(0)/l]^  ln(I^/I),  for  particular  correlation  functions 

and  values  of  , 1,1. 

0 1 o 

The  results  for  strong  scintillations  (Run  A),  a fade  level  of 
0.5  I,  and  six  pilot  levels  are  shown  in  Fig.  7.  It  is  apparent  that, 
for  this  particular  correlation  function,  indicative  of  a low  wind  con- 
dition and  strong  scintillation,  a conditioning  or  pilot  threshold  of 
0.95  I results  in  a less  than  5%  probability  that  the  irradlance  will 
fall  below  0.5  I until  18  msec  later.  If  the  conditioning  threshold  is 
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T (msec ) 

Figurt  6.  Empirical,  two-point  temporal  correlation  functions  (p)  for  log  irradiance  scintillations 


0.5  I,  the  corresponding  time  interval  is  5 msec.  As  the  threshold  is 
made  lower,  it  will  be  reached  more  often,  but  the  following  time  interval 
for  favorable  propagation  conditions  will  be  correspondingly  shorter. 

For  instance,  for  a threshold  of  0.7  I,  which  is  exceeded  7A%  of  the  time, 
and  a time  interval  of  1 msec,  the  fade  level  will  fall  below  0.5  Y only 
once  in  1,000  times.  This  approach  obviously  yields  sufficient  infor- 
mation for  an  analysis  of  use-rate  vs.  reliability. 

B.  Test  of  the  Bivariate  Normal  Distribution 


The  degree  to  which  the  empirical  scintillation  data  agree  with 
the  above  assum;>t'on  of  two-point  bivariate  normality  for  the  log  irra- 
diance  is  being  investigated  for  a number  of  runs.  For  example,  for  the 
weak  scintillation  case  of  Run  C,  the  one-point  log  distribution  (Fig.  8) 
has  an  observed  skewness  of  0.05  and  kurtosis  of  3.1.  Table  111  shows 
the  skewness  and  kurtosis  of  the  conditional  distributions  vs.  x and 


the  pilot  level  where  y and  IL  have  been  normalized  by  the  standard 


deviation  o^.  The  distributions  indicate  normality  except  for  small  x 


and  _o  = -1.  Three  of  these  conditional  distributions  are  shown  in  Fig.  9. 
a 

0 Using  measured  values  of  the  correlation  function  for  this  run 
(Table  IV),  we  compare  the  predicted,  normalized  conditional  mean  and 
variance  (Eq,  8)  witn  observed  values,  in  Tables  V and  VI.  Excellent 
agreement  with  the  gaussian  assumption  is  evident. 

Under  conditions  of  strong  scintillation,  e.g.  Run  A above,  there 
is  evidence  of  appreciable  departure  from  the  bivariate  normal  distribu- 
tion. This  is  under  further  investigation. 


IV. 


NEW  EXPERIMENTAL  WORK  AND  FACILITIES 

^ • Target  Irradiance  Statistics  with  Beam  Wander  Cancellation 
In  preceding  reports  on  this  program  (Refs.  1,5),  analytical 


results  were  given  for  the  mean  target  irradiance  and  its  fluctuations 
resulting  from  the  effects  of  turbulence  on  a finite  beam  wave,  including 
the  cancellation  of  turbulence-induced  beam  wander.  Experimental  results 
were  given  for  a 1.6  km  path  operating  at  63288  with  a wander-cancelling 
tracker . 
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TABLE  III 


> 


I 

t 


I 

Skewness  (S)  and  Kurtosis  (K)  of  Pr(x(t  + t)|x(t))  for  Run  C 


S 


\ T(ms) 

] 

2 

4 

8 

16 

0 

1 

1 

-1 

.29 

. 13 

.2 

.022 

.012 

0 

.03 

.051 

.056 

.061 

.068 

1 

-.  12 

-.017 

.01 

.13 

-.055 

2 

-.0039 

.047 

.0080 

K 

-.0065 

-.036 

\ T(ms) 

1 

2 

4 

8 

16 

0 

-1 

4.2 

3.59 

3.83 

3.29 

3.18 

0 

3.20 

3.16 

3.089 

3.05 

3.04 

1 

3.27 

3.21 

2.69 

2.77 

3.17 

2 

2.83 

2.84 

2.95 

2.98 

3.08 

( 

I 


1 


I 


I 
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I 


Conditional  I’robab i 1 i t y ("') 


- n « -1 


T * lb  mscf 


en  L(i) 


Fi^’tire  9 . Conditional  probability  for  normalized  log  irradiance 
time  t = T,  given  the  value  at  t = 0. 


* Pr  m 
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TABLE  IV 

Correlation  Function  Values  for  Run  C 

T p(t) 
2 ms  0.87 

4 ms  0.78 

8 ms  0.57 

16  ms  0.24 


TABLE  V 

Comparison  of  Predicted  and  Observed  Mean  of  Pr|x(t  + t)  x(b) 

\ o 0 

For  Run  C \ o . o 


\ 


ms 

2 4 8 16 


-1 

Predicted 

-0.87 

-0.78 

-0.57 

-0.24 

Observed 

-0.89 

-0.75 

-0.60 

-0.24 

0 

Predicted 

0 

0 

0 

0 

Observed 

-0.016 

-0.016 

-0.001 

0.008 

1 

Predicted 

0.87 

0.78 

0.57 

0.24 

Observed 

0.85 

0.72 

0 50 

0.21 

2 

Predicted 

1.74 

1.56 

1.14 

0.48 

Observed 

1.70 

1.51 

1.09 

0.45 

1 


TABLE  VI 


k 
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Comparison  of 
For  Run  C 

j\t  ms 

— \ 

Predicted 

and  Observed  Variance 

0 2 
0 

of  Prjx(t  + t) 

V ^0 

X(t)  ^0 

0 0 
0 0 

a ^ 
0 

2 

4 

8 

16 

-1  Predicted 

0.243 

0.392 

0.675 

0.942 

Observed 

0.239 

0.414 

0.651 

0.936 

0 Predicted 

0.243 

0.392 

0.675 

0.942 

Observed 

0.236 

0.406 

0.686 

0.961 

1 Predicted 

0.243 

0.392 

0.675 

0.942 

Observed 

0.278 

0.402 

0.77 

1.01 

2 Predicted 

0.243 

0.392 

0.675 

0.942 

Observed 

0.254 

0.439 

0.765 

0.99 

f1 


I 

I 

t 


t- 

f 


[ 


Following  the  above  experimental  effort,  the  facility  was  modified 

for  taking  similar  data  over  a much  shorter  path.  This  was  done  in  order 

to  achieve  smaller  ratios  of  transmltter-dlameter-to-coherence-scale  (D/p  ) 

0 

for  stronger,  reasonably  developed  turbulence,  while  maintaining  a suf- 
ficient Fresnel  number  for  the  transmitter.  This  case  scales  to  the  infra- 
red conditions  of  interest,  and  experiments  with  this  facility  have  been 
underway  for  several  months. 

The  facility  is  currently  undergoing  a further  modification,  to  add 
a simultaneous,  coincident-beam  capabilitv  at  10. 6p.  The  tracker  con- 
tinues to  operate  off  a target-point  beacon  at  6328X,  but  will  now  control 
both  optical  and  lOp  beams  (Fig.  10).  It  is  expected  that  the  turbulence- 
induced  wander,  being  a basically  geometrical-optics  phenomenon,^  will  not 
undergo  significant  wavelength  dispersion,  so  that  this  tracking  scheme 
will  effectively  remove  wander  at  the  infrared  wavelength.  Operation  at 
10. 6p  will  afford  a wavelength-scaled  test  of  our  analytical  results, 
and  will  further  reduce  D/p  in  a well  developed  turbulence.  The  receiver 
configuration  is  shown  in  Fig.  11,  and  in  the  photograph  of  Fig,  12. 

The  data  capability  and  experimental  parameters  of  this  enlarged 
system  are  summarized  in  Table  VII.  We  also  mention  that  the  analysis 
given  in  the  preceding  report,^  for  the  on-axis  mean  irradiance  in  the 
liinit  ^ ®,  is  being  extended  to  include  some  details  of  the  complete 

curve  ° \ ^ • This  is  being  accomplished  by  means  of  the 

extended  Huygens-  Fresnel  formulation  and  an  asymptotic  expansion  in 
inverse  powers  of  D/p^.  This  approach  is  free  of  the  assumptions  and 
approximations  of  a recent  Russian  treatment.^ 


B.  Long-Path  Saturation  System 

In  earlier  work  on  this  program,^  experiments  were  conducted 
over  the  longest,  lowest  path  which  can  be  reasonably  achieved  in  the 

1 


5.  J.  R.  Kerr  and  J.  R.  Dunphy,  "Propagation  of  Multiwavelength  Laser 
Radiation  Through  Atmospheric  Turbulence",  RADC-TR-74-183,  May  1974. 

6.  V.  A.  Banakh,  et  al , J.  Opt.  Soc . Am.  64,  516-518,  April  1974, 

i 

7.  J.  R.  Dunphy  and  J.  R.  Kerr,  J.  Opt.  Soc.  Am.  981-986,  August  1973. 
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Figure  11.  Diagram  of  mul t iwave length  receiver  for  beam-wander-cancellation  and  target-i 
diance  experiments. 


TABLE  VII 


Experimental  Data  and  Parameters  for  Multiwavelength  Target-Irradiance 
System 


Tracker : 
Target : 


Processing: 


Wander  signals 

Linear  irradiance  fluctuations 

Mean  linear  irradiance 

Log  irradiance  fluctuations 

and  microthermal  fluctuations 
Meterological  conditions 

Statistics  and  spectra  of  fluctuations 


Wavelength:  632SX  and  10. 6p  (simultaneous,  coincident  beams) 

Path:  91  m,  homogeneous 

Beam  height:  1.8  m 

Mean  irradiance:  120  s averaging 

Variance  1 Linear  irradiance:  120  s averaging 

\ Log  amplitude 
( Wander  angle 

Spectra:  Digital;  0.1  Hz  resolution;  100  s averaging 

Analog;  1.0  Hz  resolution;  100  s averaging 
Probability  distributions:  Digital 

Digital  sampling  rate:  1 kHz 

Target  receiver:  Quasi-Point-Aperture  at  each  wavelength 

Target  receiver  bandwidth:  1 kHz 

Target  receiver  dynamic  range:  80  db  at  63288;  60db  at  10. 6p 

Tracking  transmitter  aperture:  15.2  cm 

Tracking  transmitter  focus  resolution:  2.5ym 

Steering  servo  resonant  frequency;  300  Hz,  no  load 
Steering  servo  linearity:  <1%  of  peak-to-peak  deflection 

Steering  servo  repeatability:  0.050%  of  peak-to-peak  deflection 

Open  loop  transfer  function:  Single  Integration 

Closed  loop  bandwidth:  20  Hz 

Noise  limited  tracking  angle:  <2%  of  transmitter  diffraction-limited  width 

(63288) 

Transmitter  Chopping:  9kHz  at  10. 6y  (direct-coupled  at  63288) 

Microthermal  probe  averaging  time:  300  s 
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presence  of  earth  curvature  and  atmospheric  refraction.  In  this  effort, 
which  was  undertaken  for  a better  understanding  of  strong  turbulence 
scattering  ( hard  saturation"  of  scintillations),  it  was  determined  that 
multlkllohertz  fluctuations  and  correspondingly  small  amplitude  patches 
exist  in  the  scattered  field,  in  agreement  with  concurrent  theoretical 

g 

developments.  These  phenomena  are  also  under  intensive  study  by  Russian 
9 

investigators . 

We  have  re-established  the  long-path  facility,  and  are  currently 

completing  modifications  on  the  original  system  which  will  permit  a much 

larger  re<"elver  bandwidth  and  very  small  apertures  (<<p  ),  in  order  to 

0 

fully  accommodate  the  spatial  and  temporal  characteristics  of  the  field. 
The  system  operates  at  4880X  as  before,  but  with  100  kHz  chopping  using 
a novel  photoelastic  modulator  (Fig.  13),  and  a final  (demodulated)  band- 
width of  10  kHz.  Photomultipliers  are  utilized  to  achieve  a 100  dB  SNR 
with  the  small  receivers,  and  the  original  (dual)  covariance  receiver 
has  been  modified  for  high  resolution  at  small  separations.  A new 
demodulator  with  the  required  bandwidth  and  linear  phase  response  for 
undistorted  transients  has  been  designed  with  an80  dB  dynamic  range. 

The  system  parameters  are  summarized  in  Table  VIII.  Th'’  means 
utilized  to  achieve  high  resolution  of  the  covariance  are  showri  in 
Figure  14.  The  demodulator,  snown  in  Figure  15,  uses  a Motorola  Type 
1496  balanced  modulator  with  a Signetics  565  phase-locked-loop  to  supply 
the  necessary  constant-carrier  input;  the  latter  component  functions 
properly  over  a very  large  dynamic  range,  so  that  severe  fading  will  not 
Interrupt  the  system.  The  required  phase  or  transient  response  is 
obtained  in  a four-pole  Bessel  low-pass  filter. 

C.  Further  Plans 

Following  the  efforts  described  above,  the  wander-tracking 
system  will  be  modified  to  study  the  effects  of  enhanced  "angular 


8.  H.  T.  Yura,  J.  Opt.  Soc.  Am.  59-67,  January  1974 

M.  E.  Gracheva,  et  al,  "Similarity  Correlation j and  Their  Experi- 
mental Verification  in  the  Case  of  Strong  Intensity  Fluctuations  of 
the  Laser  Radiation",  preprint,  August  1973. 


9. 


TABLE  VIII 


Experimental  Parameters  for  Long-Path  System 


Transmitter  Characteristics 

Wavelength:  4880X  Power; 1 watt  cw 

Virtual-point-source  (Fresnel  number  = 0.002) 

Chopping  rate:  100  kHz 

Height : 7m 

Electronic  steering  for  diurnal  refraction  effects 


Receiver  Characteristics 

Height : 4m 

Diameter:  0.25mm-lmm 

Minimum  covariance  separation:  0.25mm 

Bandwidth:  10  kHz 

Dynamic  Range:  > 80  dB 

Processing : 

Variance  of  log  amplitude 

Covariance 

Statistics 

Spectra 
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scintillation"  from  a retroref lector . ^ ^ Also,  in  a follow-on 
effort,  the  scintillation  characteristics  including  covariance  structure 
will  be  determined  for  extended,  partially  coherent  (i.e.,  generalized) 
sources.  This  topic  pertains  to  turbulence  effects  on  "speckle"  and  has 
not  been  adequately  understood;  it  pertains  to  such  applications  as 
remote  sensing,  and  noise  in  adaptive  transmitters. 


10.  J.  P.  Hansen  and  S.  Madhu,  Applied  Optics  U_,  233-238,  February  1972. 

11.  R.  F.  Lutomirski,  Applied  Optics  J_4,  840-846,  April  1975. 
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RADC  is  the  principal  AFSC  organization  charged  with 
planning  and  executing  the  USAF  exploratory  and  advanced 
development  programs  for  information  sciences,  intelli- 
gence, command,  control  and  cotmunications  technology, 
products  and  services  oriented  to  the  needs  of  the  USAF. 
Primary  RADC  mission  areas  are  communications,  electro- 
magnetic guidance  and  control,  surveillance  of  ground 
and  aerospace  objects,  intelligence  data  collection  and 
handling,  information  system  technology,  and  electronic 
reliability,  maintainability  and  compatibility.  RADC 
has  mission  responsibility  as  assigned  by  AFSC  for  de- 
monstration and  acquisition  of  selected  subsystems  and 
systems  in  the  intelligence,  mapping,  charting,  command, 
control  and  communications  areas. 


